Correlation functions of one-dimensional anyonic fluids 
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A universal description of correlation functions of one-dimensional anyonic gapless systems in 
the low-momentum regime is presented. We point out a number of interesting features, including 
universal oscillating terms with frequency proportional to the statistical parameter and beating 
effects close to the fermion points. The results are applied to the one- dimensional anyonic Lieb- 
Liniger model and checked against the exact results in the impenetrable limit. 



in 

o 



c3 



I 

-a 
c 

o 
o 



> 



o 

o 

-i— > 



i 

C 

o 
o 



X 



One of the most striking results of modern physics 
is that in three spatial dimensions particles are either 
bosons or fermions, i.e. the wave-function of identical 
objects is either symmetrical or antisymmetrical under 
the permutations of space coordinates. Conversely, in 
lower dimensions generalized anyonic statistics, with in- 
termediate properties between bosons and fermions are 
possible [1]. Fractional braiding statistics in two dimen- 
sions are considered a key ingredient to describe frac- 
tional quantum Hall effect, motivating the large scientific 
effort to understand two dimensional generalized statis- 
tics. In one dimension (ID) collision is the only way to 
exchange particles and consequently statistics and inter- 
action are inextricably related. This is the main reason 
why in ID we can use bosonization and fcrmionization 
(Jordan- Wigner) to attack complicated problems that in 
higher dimensions are intractable. Excitations in ID flu- 
ids can then experience generalized statistics (not neces- 
sarily anyonic [2]). However only recently a systematic 
study of interacting ID anyonic systems started for dif- 
ferent physical motivations. First transport in quantum 
Hall fluids is localized on the edges and it is due to ID 
(chiral) anyons (see e.g. [3]). Second recent experiments 
in ultra-cold atoms succeeded in confining essentially ID 
interacting systems [4] and there is a proposal in this field 
to observe fractional statistics [5]. However a general 
treatment of correlation functions of interacting anyonic 
systems was still missing. In this letter we fill this gap 
giving a universal description of low-energy correlation 
functions in ID anyonic fluids. 

Anyonic statistics are described in terms of fields that 
at different points (x\ ^ x 2 ) satisfy 



oi**<*^)tf A ( X2 )tf A ( Xl ) , (1) 



where e(z) = —e(—z) = 1 for z > and e(0) = 0. k is 
called statistical parameter and equals for bosons and 
1 for fermions. Other values of k give rise to general any- 
onic statistics "interpolating" between the two familiar 
ones. 

In ID there is a natural way of defining an any- 
onic field in terms of a bosonic one [6]: ^^(x) = 
tf B {xy™S-»> dx ' p{x '\ where p(x) = tf A (x)V A (x) = 
^^(x)^ b(x), i.e. anyonic and bosonic density are the 



same (there is also a construction in terms of fermionic 
field [7] that will not be exploited in this letter). It is 
straightforward to show that if ^ b satisfies usual bosonic 
commutation relations, then ^ A follows Eqs. (1). In 
this letter we establish the behavior of correlation func- 
tions of anyonic gapless systems (quantum fluids) in the 
low energy regime following the original idea of Hal- 
dane [8] for standard statistics (here we closely follow 
the notations of the review [9]). The starting point is to 
split p(x) = p<(x) + p > {x), where p< refers to the slow 
parts (i.e. it consists of the fluctuations over distances 
3> po = (p(x))) and p > stands for the remaining "fast" 
modes. On a physical basis it is natural to replace in 
the definition of the anyonic field p(x) with p < (x), since 
fast fluctuations cancel in the integral. Then one can 
define the operator O(x) via d x Q(x) — irp < (x) and the 
previously introduced anyonic field becomes 

& A (x) = ^(aOe^C*) , (2) 

while the density is p{x) = p<{x)Y, m e 2mie(x) [8, 9]. In 
the bosonization (or equivalently conformal filed theory) 
language it is useful to introduce the operators <fi and 9 
via 6{x) = Q(x) - irp x and ¥ B {x) = [p^)] 1 / 2 ^ 1 ^^ 
which satisfy the commutation relation [</>(&), 0(x')] — 
iire(x — x')/2. An alternative (constructive) approach 
would have been to start from these definitions of 9 and 
(f> and show that the field in Eq. (2) satisfy anyonic statis- 
tics Eqs. (I), however we think that this way of intro- 
ducing the anyonic field is more pedagogical. 

All correlation functions in the low momentum regime 
can be obtained in terms of the so called vertex operators 



e im6(x) e in4 



(3) 



For example, in terms of these operators it is straightfor- 
ward to show that when \x\ p^ 1 the one-body density 
matrix (i.e. the Fourier transform of the momentum dis- 
tribution function) admits the expansion 

9l (x) = (^(a;)* A (0)) oc 

OO 

m— — oo 

(4) 
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that holds for all cyclic boundary conditions (the multi- 
plicative factor is determined only in terms of the high 
energy structure and it is not universal). However to 
proceed we need to fix the boundary conditions (BC). 
For simplicity, we now fix the bosonic field to satisfy 
periodic BC, obtaining (for a proof see e.g. [8, 9]) 

e -i(2ro+K)7re(z)/2 

(P0d(x)) 5 K + 2K 

(5) 

that lead to the main result of this letter 

^i{2rn-\- k,)tt pax i(2m-\- k)tt€(x) /2 

9i(x)=Po bm ~ (2m+K) 2 - • (6) 

m— — oo (P0(1(X)) 2 K +2K 

Here d(x) = Lsm(nx / L), b rn are unknown non-universal 
amplitudes that can be eventually fixed only through the 
exact/numerical computation of correlation functions in 
the specific model (as e.g. done for the XXZ spin chain 
[10] and for the Bose gas [11] using completely different 
approaches). The exponent K is instead universal and 
can be written as K = t/vj/vn. vn and vj are phe- 
nomenological parameters (with the dimensions of veloc- 
ity) introduced by Haldane to fully characterize the low- 
energy behavior of a quantum fluid [8] . They correspond 
to density and phase stiffness respectively. In practice 
they must be extracted from the solution of the model 
(when available), from numerics or, ultimately, from ex- 
perimental data. In terms of these parameters, the veloc- 
ity of the sound (i.e. the real speed of excitations in the 
interacting system) is v s = s/vJvn. For systems with 
Galilean invariance it holds vj = Vf, where vf is the 
Fermi velocity of a gas of spinless fermions with density 
po, and consequently K — vf/v s . 

Let us briefly discuss the role played by the BC. 
Fixing periodic ones for the field ^b{x) (and also for 
0(x)) implies that * A satisfy twisted BC [& A (x + L) = 
e iKvpoL ^ A (x)]. Requiring periodic BC for $>a leads to 
properly chosen twisted BC for $>b- Lack of space pre- 
vents us to show all these possibilities, but we note that 
when npoL equals an even integer the two conditions are 
equivalent. The appropriate BC for a given physical sit- 
uation can be quite specific and depend on electrical and 
magnetic properties of anyons (see e.g. the example in 
Ref. [12]). It is also easy to generalize the bosonic results 
of Ref. [9] with open BC and at finite temperature to the 
present case. Asymptotic different times correlations can 
be obtained as well. Note that, being the anyonic density 
equal to the bosonic one, the density-density correlation 
function g2(x) = (p(x)p(0)) coincides with Eq. (61) in 
Ref. [9] (obviously with K fixed by the anyonic dynamic). 
All these applications of the present method will appear 
elsewhere [13]. 

Let us now discuss our finding in the thermodynamic 
limit L — > oo. The sum Eq. (6) is dominated by the 
smallest power law in the denominator. For — 1 < k < 1 



this is always at m = 0, and for large enough systems 
we can neglect all the terms with m ^ 0. This results in 
a universal oscillating term e tK1Tp " x whose frequency is a 
direct measure of the statistical parameter. Conversely 
for fermions with k = 1 the terms at m = and m = — 1 
display the same power law (and also, by symmetry, b = 
6_i ). Consequently, for k close to 1 we expect beating 
effects due to two oscillations with very close frequency 
and practically the same amplitude. The relevance of 
this effect will be explicitely showed later on. 

Before considering specific applications, we make a last 
comment. The definitions of anyonic fields (1) have a pe- 
riodicity in k of period 2. However this is reflected in a 
periodicity of the system only if the dynamics (dictated 
by some Hamiltonian) has the same period. In our ap- 
proach this is encoded in the exponent K giving the real 
K-pcriodicity (e.g. in the Lieb-Linigcr model discussed 
below the periodicity is 4 [14]). 

The previous derivation is completely general and 
holds for any interacting ID anyonic fluid. To show the 
predictive power of this analysis we now consider two spe- 
cific examples: the free gas and the Lieb-Liniger model. 

The free anyonic gas. The momentum distribution 
function of a (more general) gas of free anyons has been 
studied in the grand-canonical ensamblc in Ref. [15]. 
Considering only the case defined by Eqs. (1), fixing 
the chemical potential in terms of the mean-density, and 
performing a Wick rotation we have 

e iKTT(p x-e(x)/2) 

gi(x) oc TTxTk > ( 7 ) 

[L sm(7rx/L)J K 

which is consistent with Eq. (6) with b m ^o = and 
K = k^ 1 as known from the dispersion relation [15]. 

The Lieb-Liniger anyonic gas. An exactly solvable in- 
teracting model is the anyonic Lieb-Liniger gas, defined 
by the Hamiltonian: 

+ c dxtf A (x)tf A (x)* A (x)*A(x), (8) 
Jo 

which describe N anyons of mass M on a ring of length 
L interacting through a local pairwise interaction of 
strength c (in what follow, we fix 2M = H = 1). The 
mean density is po = N/L. For k = the model reduces 
to the bosonic Lieb-Liniger [16], while for k = 1 just con- 
sists of free fermions. For c = oo it is equivalent to a gas 
of impenetrable anyons called Tonks- Girardeau gas [7]. 
As in the bosonic case, the physics just depends on the 
dimensionless parameter 7 = c/po and not on the two 
parameters separately. 

The model is solvable by means of Bethe Ansatz [6, 14, 
17]. The Fock space is spanned by a set of Bethe wave 
functions, each fully determined by a set of rapidities kj 
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j = 1, . . .N, solutions of the logarithmic Bethe equations 
(assuming anyonic periodic BC) [14] 

n k- - k 

kjL + ktt(N - 1) + 2 arctan 3 - 1 = 2nlj , (9) 

l=i c 

in which Ij are a set of integers (half-integers) quan- 
tum numbers for N odd (even), and d — cj cos(7tk/2) 
is an effective ID coupling of a bosonic equivalent model 
[6, 18]. If periodic BC are assumed for a correspond- 
ing bosonic field, the Bethe equations do not present the 
term ktt(N— 1) [6]. Thus, a part from an eventual shift in 
the momentum, the Bethe equations are the same as the 
bosonic Lieb-Liniger at coupling c': ID anyons with pe- 
riodic BC are equivalent to bosons with twisted BC and 
viceversa, consistently with the previous general analysis. 

For periodic BC, following [14], we observe that any in- 
teger value of k(N—1)/2 can be absorbed in the quantum 
number Ij and only changes the correspondence between 
states and quantum numbers: = Ij — [k(N — l)/2], 
where [•] stands for the integer part. For this reason 
it is worthy to define v = {k(N — l)/2}, where {•} 
stands for the non- integer part {y} = y — [y]. It is 
then trivial to calculate the speed of the sound [14] 
V s (j, k) = v s (j',0) + Aanv/L, where a = 1(0) for pe- 
riodic (twisted) BC. Consequently K(j,k) = vp/v s = 
2irpo/v s (-f' , 0) + 0(1 /L). In particular in the Tonks- 
Girardeau limit we always have K(-j = oo,k) = 1. For 
arbitrary 7, u s (7, 0) is not known analytically, but it is 
very easy to extract it numerically from the solution of 
the Bethe equations, see e.g. Ref. [9]. 

Exploiting the anyon-fermion mapping Santachiara et 
al. [19] gave an exact representation of the one-particle 
density matrix 51(2;) in the Tonks- Girardeau regime (c — > 
00) in terms of the determinant of a Toeplitz matrix for 
cases when periodic BC are equivalent to twisted BC (i.e. 
k is an integer multiple of 2/(N— 1)). These exact results 
must be reproduced by our findings Eq. (6) with K = 1. 
In Fig. 1 we compare the results of Ref. [19] with Eq. 
(6) for the same values of k plotted in Ref. [19]. For 
k = 1/10 (N = L = 21) and k = 1/2 (N = L = 61) we 
set b m ^o — and we fix bo as an overall normalization 
by imposing the correct value at x — L/2. The value 
k = 59 / 60 (N = L = 121) is peculiar because very close 
to the free-fermionic point k = 1, where two power-laws 
(m = 0,-1) have the same exponents, bo = &_i and 
all the other b m = 0. For k = 59/60 and N = 121 
the two power laws are still indistinguishable, thus to 
compare the exact results with Eq. (6) we fix bo = 6_i 
and all the other b m — 0. The beating effect resulting 
from two oscillations with almost the same amplitudes 
is evident as overall nodes (in the amplitude of the fast 
oscillations) in the real part of gi(x) at x/L — 1/4 and 
3/4, in contrast with the monotonic behavior for smaller 
k (see Fig. 1). In the three cases the agreement between 
exact and asymptotic results is perfect. 




FIG. 1: (Color online) Exact results for the anyonic Tonks- 
Girardeau gas (black lines) of Ref. [19] against the asymptotic 
large distance result (red lines) obtained here. 



Moreover in Ref. [19] the asymptotic behavior for large 
N of gi(x) was estimated using Fisher-Hartwig conjec- 
ture. The results obtained in this way are equivalent to 
the zero mode (m = 0) in Eq. (6) [20]. On the one hand 
our results give further evidence that the Fisher-Hartwig 
conjecture should be valid and is able to fix properly the 
constant bo that in our approach remains unknown. On 
the other hand, Eq. (6) shows that, even for very large 
L, close to the Fermi point the resonance effects cannot 
be caught by other asymptotic approaches (compare the 
accuracy of the plot in Fig. 1 for k — 59/60 with the anal- 
ogous one in Ref. [19]). However, hard-core anyons are 
only a specific case to which apply Eq. (6), on the same 
foot as the Lenard formula [21] for impenetrable bosons 
is only the simplest case described by the hydrodynamic 
treatment of Haldane [8] for standard statistics. 

Discussions. We provided a universal description 
of correlation functions of ID anyonic fluids in the 
low-momentum regime. Our main result is the zero- 
temperature one-particle density matrix given by Eq. 
(6). This formula unambiguously shows how to connect 
smoothly correlations of interacting bosons and fermions 
through the continuous parameter k. 

We specifically applied our method to the Lieb-Liniger 
anyonic gas, but it would be interesting to understand the 
consequences of our findings in more complicated fluids, 
as for example a generalized Tomanaga-Luttinger liquid 
[13], fractional generalization of Heisenberg spin chains 
as, e.g., the golden chain [22], fluids with internal degrees 
of freedom (e.g. anyonic extensions of the model in [23]) 
etc.. 

We find that the main universal signature of anyonic 
statistics in the correlation functions is the oscillating 
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complex term that can lead to clear signals in transport 
experiments [3]. The remaining part in Eq. (6) is very 
similar to the bosonic counterpart [8] with the general- 
ized statistics leading to a predictable "renormalization" 
of the power-law behavior for large x. This is a man- 
ifestation of the ID interplay between interaction and 
statistics in the low momentum regime. It would be ex- 
tremely interesting to understand to which extent some 
analogous correspondence can be valid for all momenta 
and in gapped systems. 
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